Observation of Exceptional Points in Electronic Circuits 
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Two damped coupled oscillators have been used to demonstrate the occurrence of exceptional 
points in a purely classical system. The implementation was achieved with electronic circuits in 
the kHz-range. The experimental results perfectly match the mathematical predictions at the ex- 
ceptional point. A discussion about the universal occurrence of exceptional points - connecting 
dissipation with spatial orientation - concludes the paper. 



PACS numbers; 03.65.Vf, 02.30.-f, 41.20.-q 
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A surprising phenomenon occurring in systems de- 
scribed by non-hermitian Hamiltonians has been ob- 
served in a number of experiments: the coalescence of 
two eigenmodes. If the system depends on some interac- 
tion parameter A , the value Aep at which the coalescence 
occurs is called an exceptional point (EP) 0. At an EP, 
the eigenvalues and eigenvectors show branch point sin- 
gularities 0, 0, S S 13 as functions of A . This stands 
in sharp contrast to two- fold degeneracies, where no sin- 
gularity but rather a diabolic point ^ occurs. EPs have 
been described in laser induced ionization of atoms 0] , in 
acoustical systems 0, and have actually been observed 
in microwave cavities 0, ITol llll. in optical properties 
of certain absorptive media [l^, uM > and in "crystals of 
light" 14] . The broad variety of physical systems show- 
ing EPs indicates that their occurrence is generic. So far 
EPs have been analyzed for the special case of a com- 
plex symmetric effective Hamiltonian |l5j used for the 
description of dissipative wave mechanical systems such 
as in microwave cavities. 

The physical interest in EPs is not only due to their 
universal occurrence in virtually all problems of matrix 
diagonalization. It is in particular the chiral character 
associated with the wave functions at the EP 16|. This 
has been experimentally confirmed recently 11] and also 
thoroughly discussed in optics for anisotropic absorptive 
media Il7| . The fascinating point is the the wave function 
at the EP: not only is there - at the point of the two 
coalescing energies - one and only one state vector, but 
- for a complex symmetric matrix - its form is given by 
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Only one EP of the pair is accessible in the laboratory as 
only one is associated with a negative sign of the imagi- 
nary part of the energy. The labels of the non- vanishing 
components in Eq. ^ are the labels of the two coalescing 
energies. Note that there is no parameter dependence; 
the form of the state vector is robust. We may associate 
this form with a circularly polarized wave. For a gen- 



eral non-hermitian matrix the ratio of the two relevant 
components is no longer ±i but an arbitrary complex 
number that can in general be associated with an elliptic 
polarisation [T^ . 

In this letter we report the results of a plain classical 
experiment where an EP has been identified and so has 
the orientation that goes with it. It is the case of two 
coupled damped oscillators. This has been proposed in 
[isl ] for two mechanical oscillators. In the quoted paper 
a discussion of the EPs is extended to general matrices 
being no longer complex symmetric. In fact, the classical 
equations of motion for the momenta pj and coordinates 
Qj read for the two mechanical oscillators 
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where ujj — ikj, j — 1,2 are essentially the damped fre- 
quencies without coupling and / and g are the coupling 
spring constant and damping of the coupling, respec- 
tively. The driving force is assumed to be oscillatory with 
one single frequency and acting on each particle with am- 
plitude Cj . Here we are interested only in the stationary 
solution being the solution of the inhomogeneous equa- 
tion which reads 

/cA 

= (iuj-M)-^ 'Jj exp(iwi). (4) 

Resonances occur for the real values u of the complex 
solutions of the secular equation 

det\iLU-M\=0 (5) 

and EPs occur for the complex values w where 
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(6) 



2 



is fulfilled together with Eq.|(51)- Note that is non- 
symmetric. 

For easy implementation electronic circuits are used 
rather than the mechanical oscillators. The momenta 
and coordinates are replaced by voltage and current, re- 
spectively. The unperturbed undamped frequencies are 
given by — l/{LjCj) with L and C being the induc- 
tance and capacitance, while damping and coupling are 
effected by resisters and mutual inductance. To allow for 
complex values of the effective inductance, resisters have 
been used in parallel to the inductance. 

The actual experimental set-up is shown in Figure ^ 
The two circuit parameters that can be changed are the 
circuit capacitance Cp and the parallel resistance Rp. 
Standard electronic components were used to construct 
the circuit. The component values are summarized in Ta- 
bleQ] The two inductors consist of multi-layer cylindrical 
air core windings wound onto PVC plastic coil formers. 
The coils are placed next to each other with their geo- 
metric axes aligned. The distance between them can be 
varied in order to change the coupling, i.e. the mutual 
inductance of the two coils. Input voltage signals va 
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FIG. 1: Measurement set-up. 



TABLE I: Component Values 



Component 


Value 


Unit 


Tolerance 


a 

Op 


65 


nP 


5% 


Ri 


3.0 


n 


10% 


Rp^ 


520 


n 


5% 


Lp 


1.80 


mH 


5% 


Ls 


2.55 


mH 


5% 


Rs 


10.0 


kQ 


5% 


R2 


3.0 


n 


10% 


Cs 


47 


nF 


5% 


M 


O.U^LpLs = 0.3 


mH 


5% 



"Can be varied from 39.0 to 77.6 riF in steps of 0.22 nF. 
''Can be varied from 430 to 610 C in steps of 10 Q. 



and vb are connected to ports A and B of the circuit 
and are used to excite the system. The circuit response 
is given by the input currents za and is into ports A 
and respectively. The two input currents can be de- 
termined by measuring the voltage drop over the series 
resistances Ri and i?2. Using Ohm's Law the currents 
can be calculated. 



The input ports A and B of the circuit have to be 
excited by ideal voltage sources. In a practical set- 
up the excitation voltage sources must have output 
impedances, Za and Zb, that are much lower than the 
input impedances of the two circuit ports. The minimum 
input impedance of ports A and B are approximately i?i 
and i?2, respectively. Buffers are used to reduce the sig- 
nal impedance, such that Za <S -Ri and Zb ^ -R2- 

A microprocessor is used to enable automatic multiple 
data collection. The microprocessor interfaces to a PC 
and varies the values of two circuit components Cp and 
Rp. The input currents are measured with a digital os- 
cilloscope. The oscilloscope is connected to the same PC 
as the microprocessor. The data from the oscilloscope 
is read directly into the PC where it is processed. The 
complete experiment is controlled with a software pack- 
age called Lab VIEW. A program written in Lab VIEW is 
responsible for communicating with the microprocessor 
and the digital oscilloscope and processing the measure- 
ment data. 

The first step in the measurement process is to de- 
termine the system Eigenvalues. The Eigenvalues can 
be determined from the system's frequency response, i.e. 
the system's spectrum. The frequency response Fa{uj) of 
the system at port A is given by ,19j 
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with AI < N. Here Va{i^) and are the spectra of 

the applied voltage signal and the input current response 
at port A, respectively. A similar expression holds for the 
frequency response at port B. For the circuit in Figure 
n it is = M = 4. The resonance frequencies are the 
zeros of the denominator in Eq.||7|). 

The system's frequency response can be determined 
by either performing measurements in the frequency or 
time domain. Time domain measurements yield more ac- 
curate results compared to frequency domain measure- 
ments. For a time domain measurement the system is 
excited by a short voltage pulse. The resulting response 
of the input current at a port is the impulse response 
of the system for that port. The Fourier transform of 
the impulse response is equal to the system's frequency 
response. The frequency response depends on the mea- 
surement configuration, i.e. the port where the impulse 
response is measured and the relative amplitude of the 
excitation pulse at ports A and B. In Eq.(I7|) only the 
roots of the numerator are dependent on the measure- 
ment configuration. The roots of the denominator, i.e. 
the eigenvalues, are independent of different measure- 
ment configurations. 

In a practical set-up it is not possible to generate in- 
finitely short pulses. However, the pulse width tp must 
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be small compared to the response time of the system. 
A valid impulse response must obey 
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with max(|a;i|) the eigenvalue with the largest modulus. 

For this experiment both ports A and B where excited 
by the same voltage pulse. The input current into port A 
was measured. Applying a Fourier transform to the mea- 
sured impulse response current the frequency response of 
the system is obtained. By fitting Eq. Q to the measured 
data the coefficients and bi can be determined and the 
resulting eigenvalues of the system calculated. 

The measured eigenvalues for different values of Rp 
and Cp are plotted on the complex energy plane in Fig- 
ure 121 Only two of the four eigenvalues of the system 
are plotted for each Rp and Cp value pair. The other 
two eigenvalues are mirrored with respect to the imagi- 
nary axis. Three resistance values were used, Rp = 430, 
Rp = 470 and Rp = 510. For each resistance value the 
capacitance was changed from 57.0 nF to 72.0 nF in steps 
of 0.22 nF. The resulting locus of eigenvalues repel from 
a point in the energy plane, which is the singular point 
of the system. The direction of each locus for an increas- 
ing value of Cp is indicated in Figure |21 We denote the 
value of the frequency at the singular point by 5*1 and 
find 5*1 « 92000 — i 11500 as indicated by a cross in Fig- 
ure 121 A second singular point 5*2 (not shown) is also 
measured and is mirrored with respect to the imaginary 
axis of the complex energy plane, i.e. iSi — (iS'2)*. The 
measurement error is less than 6%. 
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FIG. 2: Measured system eigenvalues. 

In order to measure the system phase, the system has 
to be excited with a steady-state sinusoidal signal. A 
real sinusoidal signal always consists of two frequency 
components uJi and UJ2 for which iwi = (1^2)*- Hence, 
to measure the system phase at the singular point, the 
system is excited with two complex frequencies equal to 
the two singular point values Si and 5*2. Furthermore, 
the circuit values of Rp and Cp are set to obtain eigenval- 
ues as close as possible to the two singular points. The 



excitation signals at port A and B have the form: 

v{t) = Csin(woi + 0i,)e~'^°* , (9) 

with ujo = |5R(5i)| = |5R(52)| and 70 = -3(5i) = 
—3(5*2). The phase of the signal is and the ampli- 
tude C. For the measured system the value of 70 is pos- 
itive (7 « 11500 s^^). This implies a decaying sinusoidal 
signal. 

The aim of the experiment is to measure the steady- 
state phase difference A(j)i between the two input currents 
iA and ib- This is achieved by measuring the steady- 
state current response at the two ports simultaneously. 
The system is linear and the current response will have 
the same form as the excitation voltage in Eq.©. The 
phase of the two input currents, (j)iA and (t>iB, can be 
extracted by fitting Eq.© to each of the measured cur- 
rent responses. The amplitude and phase are the fit- 
ting parameters. The phase difference is then simply 

However, the signal in Eq.© is not a power signal, 
i.e. it does not have a finite power content. It is therefore 
not possible in practice to generate such a signal. To per- 
form the phase measurement the excitation is switched 
on at t = to- For t < to the signal is zero. The on-set 
of the signal at t = to introduces a decaying transient 
response. A reliable measurement of the steady-state 
phase can only be performed at t — to 7(7^' where 
the transient response has decayed and can be neglected 
compared to the steady-state response. 

To generate signals as described by Eq.© the circuit 
values of the circuit in Figure ^ were so chosen that sin- 
gular points and eigenvalues are obtained within the au- 
dible frequency range. This enables the excitation signals 
to be generated with a standard computer sound card. 
Lab VIEW is again used to drive the sound card of the 
PC and generate any arbitrary signal within the audible 
frequency range (20 Hz to 20 kHz). The stereo output 
of a computer sound card can be used to generate two 
different signals simultaneously. 

Finally, the steady-state phase difference A0i between 
the input currents was measured while exciting both 
ports. The phase difference A4>^ between the two excita- 
tion signals was varied from to 27r radians. In Figurej^l 
the measured current phase difference A(f>i is plotted as 
a function of the excitation phase difference Acj)^. Also 
shown is the phase difference between the two input cur- 
rents, iA and ib, and the excitation voltage vb at port 
B. The measured phase difference A(j)i is largely indepen- 
dent of the excitation phase difference A(j)y with a value 
of approximately ^ and a measurement error of 5%. 

As far as the results are concerned we emphasize in 
particular two points: (i) the similarity of Fig[21with the 
corresponding figure for instance in J] or 10] is striking; 
here it refers to electronic circuits, in the quoted papers 
to atomic physics or microwave cavities, respectively; (ii) 
the measured phase difference 7r/2 which is independent 
of the phase of the excitation, underlines once again the 
universal and robust chiral character of the state vector 
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FIG. 3: Measured system phase. 



at the EP. We mention, that owing to the non-symmetric 
nature of the matrix considered (see Eq.Q) there is in 
fact a shght deviation [l8^ . For the parameters chosen the 
precise quotient of the amphtudes actuaUy is —0.08 + 
il.99 . . . and not i as in Eq.lQJ. In other words, there 
is a shght deviation from the circular toward an eUiptic 
polarization. Yet an orientation prevails. 

In this, like in all other experiments where EPs have 
been observed, dissipation is playing a foremost role. In 
fact, in quantum mechanics or optics, only absorption 
enables access in the laboratory to the singularity in the 
complex plane. And even in the experiment described in 
the present paper, where the underlying matrices are no 
longer symmetric, dissipation is crucial to give rise to an 
EP. The decay of the states provides a direction of the 



time axis. In classical and quantum systems absorption 
- and decoherence in the latter case |23| ~ occurs due 
to the interaction with the environment, the presence 
of the open channels is effectively described by complex 
interaction parameters 15. 2]| or, as in the present case, 
explicitly by Ohmic resisters. 

In turn, the wave function at the EP is usually chiral 
in character as has been discussed in |16l | and confirmed 
experimentally the first time in |ll|. We see here an in- 
trinsic connection - provided by the mathematics of an 
EP - between the direction of time and a given direction 
in space. In fact, the present experiment can serve to 
define a direction in one dimensional space. At the EP, 
it is the one oscillator whose phase is leading while the 
phase of the other is lagging. These roles are fixed un- 
equivocally for a given parameter set, unrelated to space. 
Viewing the oscillators placed along a straight line an ori- 
entation is defined along this line. This is not trivial and 
is caused by the mechanism at an EP: one, and only one- 
and-the-same oscillator has the leading phase, and this 
feature is cast in stone by a system that has no a priory 
spatial orientation. It is dissipation that brings about 
spatial orientation. 

This does not mean, at this stage, that a classical set- 
up can provide a distinction between left and right. One 
is tempted to find a truly three-dimensional setting where 
one EP might do just that. 
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